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ABSTRACT 

Context. To answer questions on the start and duration of the epoch of reionisation, periods of galaxy mergers and properties of 
other cosmological encounters, the cosmic star formation history, p, or CSFH, is of fundamental importance. Using the association of 
long-duration gamma-ray bursts (LGRBs) with the death of massive stars and their ultra-luminous nature (> lO'-erg s"'), the CSFH 
can be probed to higher redshifts than current conventional methods. Unfortunately, no consensus has been reached on the manner in 
which the LGRBs rate, pg^b or LGRBR, traces the CSFH, leaving many of the questions mentioned mostly unexplored by this method. 
Aims. Observations by the gamma-ray burst near-infrared detector (GROND) over the past 4 years have, for the first time, acquired 
highly complete LGRB sample. Driven by these completeness levels and new evidence of LGRBs also occurring in more massive and 
metal rich galaxies than previously thought, the possible biases of the Pgrb-p* connection are investigated over a large range of galaxy 
properties. 

Methods. The CSFH is modelled using empirical fits to the galaxy mass function and galaxy star formation rates. Biasing the CSFH 
by means of metallicity cuts, mass range boundaries, and other unknown redshift dependencies of the form pg^b P.(l + zf, a pg^b is 
generated and compared to the highly complete GROND LGRB sample. 

Results. It is found that there is no strong preference for a metallicity cut or fixed galaxy mass boundaries and that there are no 
unknown redshift effects (S = 0), in contrast to previous work which suggest values of Z/Zq ~ 0. 1 - 0.3. From the best-fit models 
obtained, we predict that ~ 1.2% of the LGRB burst sample exists above z = 6. 

Conclusions. The linear relationship between pg,b and p, suggested by our results implies that redshift biases present in previous 
LGRB samples significantly affect the inferred dependencies of LGRBs on their host galaxy properties. Such biases can lead to, for 
example, an interpretation of metallicity limitations and evolving LGRB luminosity functions. 

Key words. Gamma-ray burst: general - Cosmology: miscellaneous 



1. Introduction 



Long-duration gamma-ray bursts (LGRBs) are among the most 
luminous and energetic events to occur in our Universe. First 



signalled by their prompt high ene rgy emission (observer frame 
y-rays; e.g., iNaravan et al.l 11992 ). pinpointing regions of star 
formation irrespective of host galaxy luminosity, they are then 
followed by a longer wavele ngth and longer lasting afterglow 
(X-ray through to radio; e.g., Cavallo & ReesI |1978; Goodrnag 
ll986HPaczvnskilll986l: lM eszaro sl l2002h from which precise po- 
sitions, redshifts and host galaxy abundance measurements can 
be obtained. Such high luminosities allow LGRBs to be de- 
tected to high redshifts making them powerful probes of the 
early Universe. As early as 1993, they were thought to have 
been the result of the core-collapse of a mas sive star, under 
the mechanism named the collapsar model (WooslevI [1991 
[Paczvnski 1998; MacFadyen & Wooslev 1999). Given strong 
evidence of an association o f a LGRB to a supernova (SN) 
in 1998, (lGalamaetal.lll998l SN1998bw) and first conclusive 
spectros copic confirrnation of SN2003dh associat ed with GRB 
030329 (IStaneket al.ll2003l: iMatheson et ani2003h . the connec- 
tion of most LGRBs with the core-collapse of massive stars is 
now unquestionable (for a review, see e.g., IWooslev & Bloonil 
I2006h . 



As a result of the collapsar model's connection to the death 
of massive stars and the small distances that the progenitors 
travel away from their birth location in their relatively short 
lives, LGRBs can be used to trace the star formation rate of 
their host environment. Combined with the fact that LGRBs 
have been spectroscopically confirmed to be a cosmological 
phenomena (see e.g., Metzgeret al. 1997), they could be used 
to trace the cosmic star formation history (CSFH; p.). The 
advantage of using LGRBs rather than conventional methods 
based on Ly man-break galaxies ( LBGs) and Lyman-a emitters 
(LAEs) (e.g.. lRafelski"etani20Tl ). is that their immense lumi- 
nosities allo w them to be detecte d up to very high reds hifts, 
e.g.. z ^ 6.7 dGreiner et alJl2009ah . 8.2 dTanvir et al.ll2009h and 
9.2 (ICucchiaraetal .1120111) . Therefore, the CSFH-LGRB rate 
(LGRBR) con nection has been investigated from as early as 
1997 (see e.g.. lTotanilll997HWiiers et alll998h . 

Knowing the CSFH to such high redshifts is of fundamental 
importance for studying galaxy evolution. The current picture 
of the C SFH ( s ee e.g., Madau et al.lll996t iHopkins & BeacomI 
I2006t 12001 iKistler et al.. ,2009ir is that there is a steady in- 
crease of star formation from z = to z = 1 followed by a 
plateau up to redshift z ~ 4. However, the shape of the CSFH 
above z = 4 remains highly uncertain, where it could continue to 
plateau fe istler et al l2009h . drop off (lLil l2008h or even mcrease 
(iDaigne et al.ii2006r 
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The achievement of constraining the CSFH to high red- 
shifts would allow many questions about the Universe to be 
answered, with one of the most sought after being: "when did 
reionisation occur?". To study such questions at high redshift 
an understanding of the relationship between the CSFH and 
LGRBR would be incredibly valuable. However one obstacle 
that plagues the CSFH-LGRB rate connection is whether or not 
LGRBs are biased tracers of star formation. The question of bi- 
asing was first introduced because core-collapse models could 
not generate a LGRB witho ut the prog enitor sys tem having low 



signal-to-noise). This has facilitated, for the first time, highly 
complete GRB redshift samples. 

Only recently have highly complete samples been used to 
show that GRBs that exhibit no afterglow are primarily the re- 
sult of t hese GRBs o riginating in a galaxy of hig h extinction 
( Greiner et al. 2011a ). These dark bursts (see also iFvnbo et al] 



2009a; Cenko et al.1 120091: iPerlev et al.l |2009') are usually not 



included in GRB host follow-up programs due to poor lo- 
calisations. New evidence has shown that GRBs with heav- 
ily dust-extinguished afterglows exist in systematically more 



metal l icity 0.3 Z^; e g.. iHirschi et al. 2005: lYoon & Laiige^ massive galaxies (Kruhler et al. 2011a), thus originally bias 



l2005t IWooslev&Hegeii 12006V Also, it was noticed that the 
LGRBR was flatter at higher redshifts than the CSFH and 
that LGRBs were t ypically found in low - mass, low-metalLicity 
galaxies (see e.g., Le Floc'h et al.' '2 0031: Savaglio et al.l l2009t: 
ISvensson et al.ll2010t Mannucci et al.l 1201 1 ). The difference in 
the LGRBR and CSFH, suggested in previous papers, can be 
understood if cosmic metallicity thresholds, evolving LGRB lu- 
minos ity distributions (e .g., the progenitor evolves with red- 
shift; Virgili et al. 20n), evolving stellar initial mass func- 
tions^ (Wang & Dai 201 1 ) o r samp le selection effects (see e.g., 
ICoward et aLlfeOOSHLiTet al.ll2012l) are included. 

A significant hindrance of the analysis of the CSFH-LGRBR 
connection is the number of LGRBs with no redshift mea- 
surement, i.e., redshift incompleteness. For a redshift to be 
measured, high precision localisations (~0.5") are required for 
foflow-up spectrosc opy. Since the advent of the Swift satellite 
dGehrels et al.ll200 4l). a survey telescope equipped with a GRB 
alert telescope (BAT), an X-ray telescope (XRT) and a ultra- 
violet/optical telescope (UVOT), there have been over 600 GRB 
detections in the past 7 years. The combination of large sky 
coverage of the BAT and the (near-) instantaneous follow-up 
with the XRT and UVOT, high precision localisations have been 
achievable, which were not as common in the preceding era. 
Accompanying this, the number of robotic ground-based follow- 
up telescopes has increased , allowing GRBs to be seen minutes 
and even seconds (Rykoff et al.l2009l) after they trigger the BAT. 
However, despite such efforts it has only been possible to obtain 
redshifts for ~ 30% of GRBs in comparison to the high success 
(84%) of X-ray detection^ 

Such a low redshift completeness is the result of: large un- 
certainties in the GRB localisations, weather, the GRB does not 
fit the redshift follow-up program criteria and GRB sky location, 
to name but a few. This biasing can become so complex, that 
it is not always so simple to remove from the s ample consid- 
ered , but none the less has been tried before (e.g.. lCoward et all 
12008 ). Consequently, many groups have been trying to im- 
prove the completeness levels of their statistics. One such in- 
strument has been set to this task, the gamma-ray burst opti- 
cal near-infrared (NIR) detector (GROND; Greiner et al. 2008). 
This multi-channel imager mounted at the 2.2m MPG/ESO tele- 
scope at La Silla (Chile) has operated for the past 4 years as 
an automated GRB afterglow follow-up instrument. The robotic 
nature of GROND has allowed it to significantly increase the af- 
terglow detection rate (~ 98%) for all LGRBs observed within 
4 hours of the trigger Due to the reduced attenuation from gas 
and dust at observer frame NIR wavelengths, GROND's multi- 
band capabilities in combination with its rapid-response allows 
photometric redshifts to be determined when spectroscopic ob- 
servations were not possible or fruitless (i.e., redshift desert, low 



ing our opinion on host galaxies. Furthermore, it has been 
seen before that not all GRB host g alaxies fit the picture 



Levesque et al 



picture 

T l2oTot 



of low-mass, low-metalLicity (see e.g.. .^•^y^^^^^ ^.^.^ .^^^^i. 
[Hashimoto et al.ll2010l: ISava glio et al. 2012^) and that some can 
extremely dust extinguished galaxies (L evan et ^ 



in extremely dust extinguished galaxies (Lei 
iBerger et all 120071: [Hashimoto et al.l 120101: |H 



unt et al 



occur 
2006; 

201 U Rossi et al. 11), which are also usually associated with 
large metallicities. 

In this work we investigate, simultaneously, the effect of 
mass ranges and metallicity cuts placed on the CSFH and 
its effects on the CSFH-LGRBR connection wit h the recently 
availab le and highly redshift-complete sample of iGreiner et al.l 
(1201 lah . The paper is outlined as follows. The CSFH and 
LGRBR models are explained in Sect. |2] The description of in- 
dividual LGRBs and their properties are discussed in Sect. [3] 
The implementation of the model with the data is presented in 
Sect.m The LGRB host property results and differences to other 
studies are discussed in Sect.|5] In Sect. |6] we discuss the caveats 
of our approach and make predictions on the LGRB distribution 
and CSFH at high redshifts, and summarise our conclusions in 
Sect. |2] A ACDM cosmology has been assumed throughout this 
paper (Qa = 0.7, Qm = 0.3, Hq = 73.0 km s"' Mpc"'). 



2. CSFH & LGRB Rate Models 

The aim of this work is to compare LGRB number densities 
generated by models, that vary with regards to their host galaxy 
characteristics, to an observed LGRB number density. To do this, 
the CSFH is converted to a LGRBR utilising a conversion fac- 
tor, ?7grb (i.e., from the stars in a given galaxy how many will 
produce a LGRB that is then observed by a given instrument). 
The number density is then calculated by numerical integration. 
The whole process is outlined in the next section and is divided 
in the following way: 

A. CSFH Model 

1. The CSFH is formulated from empirically constrained 



models of galaxy star formation rates (Sect. 12.1b and 
galaxy mass functions (Sect. l272l i. 
2. Restrictions are implemented in the CSFH 
model (Sect. 12.51 1 on such things as: galaxy mass 
ranges (Sect. 12.2b . galaxy metallicities (Sect. 12.3b and 
red-dead galaxies (Sect. lZ4l i. 

LGRBR Model 



' Taken from the Swift GRB tables; 
|http : //heasarc . nasa . gov/docs/swif t/archive/grb_table/stats/ | 



1 . An initial mass function for stars is chosen (Sect. 12.6b . 
which will be used for T/^ib. 

2. A LGRB luminosity function is chosen to model the 
samples that will be investigated and implemented in 
rjgrb (Sect. 12.7b. 

3. The LGRB number density is calculated using numerical 
integration and ri^^b (Sect. lZSl l. 
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function parameters taken from 



2.1. Star Formation Rate 

During a galaxy's evolution, its gas supply will go through two 
distinct phases named the source and sink. The source phase is 
the inflow (accretion) of cold gas due to the halo's potential well, 
and the sink phase is the consumption/loss of gas by production 
of stars and outflows. These two processes occur constantly dur- 
ing the evolution of a galaxy with different weighting, e.g., at 
redshift 2 it is still possible to have cold accretion on to massive 
galaxies (see e.g., Khochfar & Silk 2009b). However, it has been 
shown via simulations that galaxies wil l reach a steady state be- 
tween these two phases (B ouche et al.ll201Q) and that the SFR 
can be considered to be solely dependent on the inflow of cold 
gas. It is this process that is believed to result in the SFR "se- 
quence", i.e., the relationship between the star formation rate 
and stellar mass, M,, of a g alaxy. This sequen ce is seen to obey 
the following relationship (iBouche et alJl2010[) : 



SFR{M„z) = 150(m./10"Mo)° 



0.8 / 1 + Z 



3.2 



2.7 



Mgyi-- 



(1) 



The main drawbacks of this formulation are that it is only fit 
up to redshifts of z = 2 (for a model that can account for 
high redshift observati ons of the SFR stellar mass relation, see 
lKliochfar&S"iiai2011h and that it is built from average stellar 
masses. This results in a spread of ~ 0.3 dex. The final model is 
compared to data in Sect. 12.51 



2.2. Galaxy Mass Function 

The galaxy mass distribution functio n (GMF) is c ommonly de- 
scribed by a Schechter function (Sch echteij|1976l) . However, as 
galaxies evolve they accrete more gas and create more stars, thus 
modifying the number density of ga laxies at a given mass for 
a specific redshift (e.g., iKhochfai^ Silk 2009a). Therefore, a 
redshift evolving GMF, as measured from the GOODS-MUSIC 
field (Fontana et al. 2006), is used in this paper (similarly to 
lYoung & Frverll2007l:iBelczvnski et al.ll2010h : 



'(M,z) = 0*(z)ln(lO)[lO^-«'('^Y""**-^'^^-'°""'", 



(2) 



where Ai = logjo {M^/Mq), M, is the stellar mass of the galaxy 
and the parametric functions obey: 



Table 1: Galaxy mass 
iFontana et aLl(l2006l) . 





m; 


M* 






11.16 


0.17±0.05 


-O.OVh 


tO.Ol -1.18 


a; 


'P'n 


^; 






-0.082±0.033 


0.0035 


-2.20H 


tO.18 



fit relation, c, is ass umed throughout this paper, taken from 
ISavaglio et all (l2005h . of flie fornS 



e[0/H] = e(M,,tH) = -7.5903 -H 2.53 15 logM. 



-0.09649 log^ M. -I- 5.1733 logfn 
-0.3944 log2 tH - 0.0403 log tn log 



(3) 



where M, is the stellar mass and tn the Hubble time measured 
in 10^ years. The relation is taken to be true over the redshift 
range we investigate (up to z ~ 12), but has been note d that it is 
difficult to reconcile for z > 3.5 dMaiohno et alj|2008l) 



2.4. Galaxy Downsizing 

Massive galaxies in the redshift regime of z ~ 1-3 have 
been seen to have little (passively evolving) or ii o (dead) star 



been seen to have little (passively evolving) or ii o (dead) star 
formatio n dCimatti et alj 120041: iLabbe et al.l 120051: iDaddi et al.l 
2005; K ong etal.ll2006^ Observationally, more massive galax- 
ies are seen to cease star formation first and this has been 
coined downsizing. The mechanism to explain this process 
is believed to be caused by the quenc hing of star forma - 
tion by active galactic nuclei (AGN; e.g.. ICroton et alj|2006l) . 
galaxy mer gers (e.g.,[S pringel et al. 2005), temperature thresh- 
olds (e.g. , Dekel & Birnboim 2006) or gravitational heating 
(e.g., iKhochfar & Ostriker 2008), but the relative importance 
is still unclear. It is possible to experimentally constrain the 
quenching stellar mass, Mq, above which no star formation is 
occurring. This is als o seen to evolve st rongly with redshift in 
the following manner dBundv et al.ll200d) : 



b* (z) = 0*,(l+z)' 



Mq(z) = M«(1+z> 



3.5 



(4) 



a* (z) - aQ + a^z 

M* (z) = M; + m;z + M^z^ 

The parameter values are given in Table [1] Given that the above 
GMF is modelled on flux limited surveys, we note that methods 
have been used to calculate missing galaxies by mass-to-light 
ratios. Secondly, the empirical fits are al so only applicable for 
redshifts, z < 4. See lFontana et all (l2006l) for fufl details on the 
galaxy selection and GMF fitting. 

2.3. Mass-Metallicity Relation 

As galaxies evolve and stars form, they pollute the galaxy 
with metals through stellar winds, supernovae and other forms 
of feedback and therefore can mediate the inflow and col- 
lapse of gas. A mass-redshift dependent metallicity empirically 



where Mq is the local (z ~ 0) quenching mass. A value of 

logio (mJ^/Mq) = 10.9 ± 0.1 dVillar et al.|[2011h is assumed 
throughout. 

2.5. CSFH 

Given the previous models, it is possible to model the CSFH. 
This is simply the integration of the individual galaxy SFRs over 
the GMF. As our aim was to also investigate varying character- 
istics of galaxies that contribute to the CSFH, we include the 
following free parameters: metallicity upper limit, and stellar 
mass upper and lower limits. Mi , M2. This results in a CSFH of 
the form, 

- The conversions 6 [0/H] + 12 = 6 [Z/Zq] + 8.69, jSavaglioll2006h 

and tn = — — 7= — are used throughout this paper. 

i.4ynA«o 
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Mode 

Hopkins & Beacom (2006) 
Coleetal. {2001} 
Hopkins et ai. (2004} 
Li {2008} 




Redshift (z) 



Fig. 1: The CSFH given by the model in Eqn.|5]using no metal- 
licity limit (i.e., no el) and (M1M2) - (7, 12) (bold-red line), 
and two empirical fits (green-d otted and blue-d ashed lines). The 
bla ck circ les denote data from lHopkinsI ( |2004|) and purple from 
lLil d20Q8h . 



((z)y(M,,z,eOSFR(M,,z)4>(M,,z) dM, 

Jm, 



(5) 



where the metallicity and downsizing constrains y and f respec- 
tively, are defined by the following relations: 



1 ife(M„z) < fL 
if e(M.,z) > fL 



ifMQ(z) > M, 
ifMqiz) < M, 



(6) 



(7) 



however, simulations we ran using no quenching mass gave the 
same results that are outlined in Sect.|5] 

The important difference between the approach used in this 
paper to others is that the form of the model in Eqn. |5] allows 
freedom for the mass and metallicity of the contributing host 
galaxies. Also, as we are interested in the low-redshift range 
(z < 3), we can remove the freedom of the parameters of the 
individual models incorporated (e.g., Eqn.[T]). 

2.6. Initial Mass Function 

Previously, it was mentioned that different types of initial mass 
functions (IMFs) have been used to explain the current LGRBR- 
CSFH connection. Cuii'ently there is no co nsensus on the correc t 
IMF to be used: Salpeter (Salpeter"l955), Scalo (Scalo" l986l) . 
Kroupa(Kroupa 2001) or Ch abrier ( Chabrier 2003), and is still a 
lively debated issue (see Ba stian et alj 2010). A commonly cho- 
sen IMF is the Salpeter IMF and is used throughout this work: 



i// (m) — m 



(8) 



where m is the star mass and a - 2.35. We note that the collap- 
sar model requires the formation of a black hole and this only 
applies to the high end mass of stars (> 3OM0). The main differ- 
ence between possible IMFs is at the low-mass end, and there- 
fore any IMF chosen would have been viable as long as it was 
redshift indepen dent. We do not consider an evolv ing IMF (top 
heavy; see e.g.. i Dav^ 120081: IWeidner et al.l 1201 Ih, but similar 
studie s have been done utilising this type (e.g., I Wang & Pail 
1201 Ih . 

2. 7. GRB Luminosity Function 

The LGRB luminosity function (LF) can take on many forms, 
for example: Schechter functions, broken power-laws, log- 
normal functions and normal functions (Gaussian function). 
Throughout this paper a normal function ha s been used (see e.g., 
iBastian et al]|2010l: iBelczvnski etal]|2010h of the form: 



' (L) - UtC 



,(L-L,)/2o-^ 



(9) 



The form of the above two step-functions (Eqns. |6] and [TT l 
could be different. For example, IWolf & Podsiadlowskil (l2007h 
implemented a metallicity cut-off of not only step-functions, but 
power laws and broken power laws. Their results show that they 
cannot currently discriminate between the form the metallicity 
cut can take. Secondly, it has been seen observationally that the 
quenching mass is not a single cut-off but a smooth function of 
mass (e.g., Kong et al. 2006), however, the knowledge of this 
trend is also not well known. Therefore, for simplicity both are 
chosen to be implemented as step-functions. 

The numerically in tegrated model, prese nted i n Fig. [1] is 
compared to data from Hopkins ( 2004) and Li (2008), and to two 
empiri cal fits from Hopkins & Beacom (2006) and ICole et alJ 
(l200l . The behaviour of the three curves is similar, peaking 
between redshift 2 and 3 and becoming negligible above red- 
shift 7. A goodness of fit of our model to the data yields a 
lA.o.f. =36/50. There is a small overproduction of stars at 
low redshift in comparison to the data, which can be lowered 
by reducing the quenching mass, Mg. However, there is a range 
of possible quenching masses that could be used and still give a 
good fit to the experimental data. As a result, i t was o pted to use 
the local quenching mass given bv lVillaretal] ( l20Ilh . It is noted 
that this parameter could also have been left free in our analysis. 



where L - log ^L,so/erg s One of the fits to our data can be 
seen in Fig.|2] and the best-fits for all the samples considered (see 
Sect. |3]l can be found in Table |2] The total fraction of LGRBs 
seen by a specific instrument, then takes the form: 



fiz) 



r 

JLr, 



>(L)dL 



(10) 



The luminosity limit, Liimit(z), of the sample can be calculated us- 
ing the luminosity distance, Dl, of the form Lnmit = ^-^D^Fw^nii- 
By taking Fiimit to be the lowest luminosity of the sample (see 
Sect. |3]l, results in the following flux limit: Fiimit = 1.08 ■ 
10"^ erg s"' cm"^. This flux limit is essentially the limiting flux 
of the BAT on board the Swift satellite. 

It is also possible that the LGRB LF is a function of 
redshift. A redshift dependence would imply that there is a 
change in the LGRB explosion mechanism throughout redshift 
or that the progenitor mass distribution, or mode (single, bi- 
nary), is changing. This has been investigated previously by 
SalvateiTa & Chincarini (20()7|) and ICampisi et al. (2010) who 
both come to the conclusion that a non-evolving LF is possible, 
but requires a metallicity cut to reproduce the observed LGRB 
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Table 2: GRB Luminosity Function best-fit parameters. 



Sample 


log,o(LJergs ') 


0"L 


11, 


x'/d.o.f. 




S 


51.33 


0.97 


10.89 


0.63/4 


P 


51.74 


0.93 


13.19 


5.56/5 


ULl 


51.06 


0.87 


23.48 


0.52/3 


UL2 


51.11 


0.92 


22.95 


0.39/3 



dist ribution (for metalli city cuts see Sect.|4]i. On the contrary to 
this, iButler et"al] (12010 ) rule out an evolving luminosity func- 
tion to the 5cr level using Swift data. To increase the simplicity 
of the modelling, we assumed that the mechanism was redshift 
independent and that the progenitor was of a constant type. For 
a discussion on implications of these assumptions, see Sect. 16.31 



2.8. Long Gamma-Ray Burst Rate 

Given the CSFH model formulated in the previous section, we 
now consider which population of galaxies contribute to the ob- 
served LGRBR. To do this, six parameters have to be considered; 
(i) the co-moving volume; (ii) the instruments limiting flux, that 
affects how many LGRBs of a given luminosity can actually be 
detected; (iii) a LGRB probability, that transforms the number of 
stars being formed to the number of LGRBs being produced; (iv) 
proportionality relation of the CSFH-LGRB rate, to include any 
further dependencies unaccounted for in our model; (v) the frac- 
tion of the sky observable by the instrument; (vi) and the length 
of time the instrument has been running. The total detectable 
LGRB rate for a given instrument in a redshift bin {z\ , zi) is then 
simply: 



Nizuzi) = ?7gi-b 



X 



■^2 f(z)p(z,eL,MuM2)il+zff 



dz, (11) 



Dcom is the co-moving distance. 



where f =47rD2„,„,., 
T/grb is the probability of stars resulting in a LGRB and then 
detected by an instrument, and 6 is the power of proportional- 
ity for t he CSFH t o LGR BR (see e.g., Bromm & Loeb 2006; 
Langer & Normanl 20061; lO aigne et al. 2006; Young & Fryer 
20071; JWolf & Podsiadlowski ' 2007; Salvaterra & Chi ncarinii 



2007t iK isder et al. 2009; Campisi et al. 2010; Oin et al.l l201( 
Butler et al. 2010; Wanderman & Piran 2010; Belczynski et al.' 
2010: Virgin etal 2011; de Souza et al. 2011; Wang & Dai 



201 it llshida et alJl201 ih . The LG RB probability can be param- 
eterised as (explained similarly in lBromm & Loebll2006h : 



Hgrh - AQ ?7coll rjBH '7tirae 'yX-iay '7redshift '7other, 



(12) 



where AT" is the length of observations, AQ. is the solid angle of 
the observable sky and the remaining parameters are discussed 
in the following paragraphs. 

As mentioned, LGRBs are likely associated with black holes 
produced by the collapse of stars above a specific progenitor 
mass, Mbh, and is quantified in the following way: 



L. °™ (m) dm 

_ JMbh ^ ^ ' 

'?BH - , 

mi^(m)dm 



(13) 



where M,nin and M^ax are the minimum and maximum star 
masses considered in our model respectively. Secondly, GRBs 
are believed to be jets of coUimated matter, as a result of breaks 



observed in afterglow light curves. This means that only a small 
fraction of bursts are visible to the observer, quantified by t^coU- 

The last four probabilities, 77time, '/x-ray, '7i-ed.shift and ?7u„known 
are a result of our sample selection criteria, outlined in the next 
section. They are the probability of detecting an optical-NIR af- 
terglow less than 4 hours after a GRB trigger, the probability of 
detecting an X-ray afterglow, the probability of detecting a red- 
shift, and any unknown probabilities respectively. 

Finally, the parameter 6 includes further dependencies be- 
tween the CSFH and LGRBR not a ccounted for by our model. 
This parameter is the same utilised in lKistler et alJ (l2009l) . where 
it is used to incorporate effects that are not known, i.e., a black- 
box approach. We employ the same idea, so that a non-zero 6 
would imply we are missing an effect in our CSFH modelling. 
One example is an evolving LGRB LF or an evolving stellar 
IMF. It is now possible to formulate a LGRB number density dis- 
tribution for a redshift bin through the use of a modelled CSFH, 
which can then be compared to an experimental data set. We 
would like to note that the empirically calibrated relations used 
in the model are only verified for redshifts z < 4, and could differ 
at higher redshifts. To investigate the implications of this on our 
final results, we manually modify the CSFH at redshifts z > 3 in 
Sect. I 



3. Data Samples 

3.1. Gamma-ray Burst Sample 

Our LGRB sample is taken from iGreiner et al! which 
unlike previous studies, is highly complete in terms of op- 
tical/NIR afterglow detection rates and measured redshifts. 
iGreiner et al.l (1201 lal) chose the sample by selecting GRBs that 
have been detected by GROND within 4 hours after the Swift 
BAT trigger and that exhibited an X-ray afterglow. This selection 
results in a sample of 43 GRBs: 39 LGRBs and 4 short GRBs 
(believed to be associated with th e merger of 2 neutron st ars or 
a neutron star-black hole svstem: feelczynski et al.ll2006l) . The 
39 LGRB sample contains 3 1 spectroscopic redshifts, 6 photo- 
metric redshift measurements (3 of which are upper limits) and 
2 with no optical/NIR afterglow detections and thus no redshift 
measurements . For more details on t he burst sample and individ- 
ual bursts, see IGreiner et al.l ( 1201 lal) . 

As mentioned previously, it is only LGRBs that are believed 
to trace the death of massive stars, and so the short-GRBs were 
not considered in our analysis. The LGRB sample was then sub- 
divided into; spectroscopic (S): bursts with spectroscopic red- 
shift (31/39); photometric (P): spectroscopic sample including 
photometric redshifts (34/39); upper limit 1 (ULl): photometric 
sample including upper limit photometric redshifts (36/39); and 
upper limit 2 (UL2): upper limit sample including a possible red- 
shift measurement (37/39). This subdivision was introduced for 
two reasons. Firstly, to see if the results changed when having 
a single method (spectroscopic) of redshift identification com- 
pared to multiple types of redshift identification (spectroscopic 
and photometric). Secondly, the upper limits were included to 
get as close to a 100% redshift-detected-sample as possible. All 
the samples cover a range of z = to z = 6.7, and the P, ULl 
and UL2 samples cover z = to z = 9.2. Any bias in the redshift 
distribution, i.e., possible deficits in the high redshift regime and 
effects of selecting bursts followed by GROND < 4 hours af - 
ter the trigger, has been investigated bv Grein er et alJ (1201 lah . 
This was carried out by comparing the sample to one of twice 
the size (Fynbo et al. 2009a), with a completeness at the time 
of 50%, utilising a Kolmogorov-Smirnov test. This showed that 
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Fig. 2: Number density of the luminosities for all of the sub- 
samples considered. The bold line depicts an example best-fit 
normal-Gaussian to the spectroscopic sample. Bin sizes are cho- 
sen for presentation of the data. 
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Fig. 3: Luminosity-redshift space of all the subsamples consid 
ered. The dashed line depicts the flux limit, found by fitting the 
spectroscopic distribution, Fu^nit - 1.08 
SectHjl. 



10 erg s ' cm ^ (see 



the difference at high-z was statistically insignificant and at the 
Icr level, both samples could be drawn from the same underly- 
ing distribution. For information on individual bursts see Table[3] 
and for the property-distributions see Figs.|2][3]and|4] 



3.2. LGRB Properties 

For each burst the measurement of a redshift and luminosity is 
required for the modelling of the CSFH-LGRBR, which will be 
discussed in Sect. |2] Redshifts fo r all of the samples w ere ob- 
tained from the table compiled bv lGreiner et al.l (1201 lah and the 
luminosities were calculated from the isotropic energy (Siso) and 
the duration of the burst over which 90% of the flux is released 

i Tqo), using the st andard relation of Liso = £i,so(l + z) /Tgo 
Bloom et al.ll200lb . The iiiso and T^q for each burst was ob- 
tained from fits to the prompt emission taken from the exten- 



Fig. 4: Number density of the redshifts for all of the subsamples 
considered. Bin sizes are chosen for presentation of the data. 



sion of lButler etall (l2007l from here NBllfl For bursts with 
photometric or upper limit reds hifts an E\^n was ca lculated via 
the same procedure described in'Bloo m et al.l (1200 lb . One of the 
following energy spectra, (E), was used based on the best-fit to 
the LGRB's BAT spectra: 



iff 



'^peak 



Power Law (PL) 
Cut-off Power Law (CPL) 



(14) 



(see e.g. ISakamoto et aT]|2008h where Ka is the normalisation 
at £■ = AkeV in units of photons keV"' cm^^s"', fipeak is the 
peak energy in the cut-off power-law spectrum, i.e., in E^<p {E) 
space, and a is the spectral slope. Utilising the best-fit model, 
the K-corrected Eiso is calculated vi£0: 



'1 — 



E4>{E) dE, 



(15) 



where Di is the luminosity distance, z the redshift, and Af the 
time over which the spectra is fit (a time-integrated spectra). 

4. Methodology 

We moved through a 3D parameter space, consisting of mass 
ranges, metallicity upper limits and proportionality relations 
(i.e., 6), in small step sizes and at each step compared the mod- 
elled LGRBR produced to our sample, using least statistics to 
asses the goodness of fit. The parameter details and methodology 
used are explained in the next sections. 

4. 1 . Fixed Parameters 

The parameters AT", Af2 and T/git in Eqn. [TT| were assumed to 
be independent of redshift and luminosity. The values used can 
be found in Table |4] but are discussed briefly in the following 
paragraphs. 



^ http : // astro . berkeley . edu/ ~ nat/swi f t/| 

For BAT it is assumed £1 = 15keV and & = 150 keV, but we 
rescale similar to NB 11 to the range of, Ei 



1 keV to £7 = 10" keV 
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Table 3: The properties of the GROND Dark Burst samples. 



Burst 


z" 


log,o(£iso/erg) 








Prompt Model 


Sample ^ 










Name 


a 












(s) 






(keV) phkeV"' cm^'^s"' (s) 




GRB100316B 


1.180 (1) 


51.57 


4.3 








S, P, ULl, UL2 


GRB091127 


0.490 (2) 


52.36 


9.6 








S, P, ULl, UL2 


GRB091029 


2.752 (3) 


52.84 


40.0 








S, P, ULl, UL2 


GRB091018 


0.971 (4) 


51.88 


4.4 








S, P, ULl, UL2 


GRB090926B 


1.24 (5) 


52.52 


126.4 








S, P, ULl, UL2 


GRB090814A 


0.696 (6) 


51.30 


113.2 








S, P, ULl, UL2 


GRB090812 


2.452 (7) 


53.89 


99.8 








S, P ULl, UL2 


GRB090519 


3.85 (8) 


53.76 


81.8 








S, P, ULl, UL2 


GRB090313 


3.375 (9) 


53.11 


90.2 








S, P, ULl, UL2 


GRB090102 


1.547 (10) 


53.65 


30.7 






... 


S, P ULl, UL2 


GRB081222 


2.77 (11) 


53.17 


33.5 








S, P ULl, UL2 


GRB081121 


2.512 (12) 


53.08 


19.4 








S, P, ULl, UL2 


GRB081029 


3.848 (13) 


53.48 


169.1 








S, P, ULl, UL2 


GRB081008 


1.968 (14) 


52.83 


199.3 






... 


S, P, ULl, UL2 


GRB081007 


0.529 (15) 


50.89 


5.6 






... 


S, P ULl, UL2 


GRB080913 


6.7 (16) 


52.72 


8.2 








S, P, ULl, UL2 


GRB080805 


1.505 (17, 18) 


52.99 


111.8 








S, P, ULl, UL2 


GRB080804 


2.20 (19, 18) 


53.94 


61.7 








S, P, ULl, UL2 


GRB080710 


0.845 (20, 18) 


52.69 


139.1 








S, P ULl, UL2 


GRB080707 


1.23 (21, 18) 


52.07 


30.3 








S, P, ULl, UL2 


GRB080605 


1.640 (22, 18) 


53.31 


19.6 








S, P, ULl, UL2 


GRB080520 


1.545 (23, 18) 


52.68 


3.0 








S, P, ULl, UL2 


GRB080413B 


1.1 (24, 18) 


52.17 


7.0 








S, P ULl, UL2 






53.41 


46. / 








C D T TT 1 T TT T 


GRB08041 1 


1.03 (26, 18) 


53.60 


58.3 








fl n T TT 1 T TT 

S, P, ULl, UL2 


GRBU8Uj3U 


1.51 (27) 


52.13 


66.1 






... 


C Tl T TT t T TT n 

o, P, ULl, ULz 


VJivDUOUZ iu 


9 f^A^ ^98 1 8^1 


J J.UD 










^ P ITT 1 ITT 


VJlvDUoU i Ay 


^.jH-y \Z.y ) 












*^ P ITT 1 ITT 9 


GRB071031 


2 692 (301 


52.91 


187.2 








S, P, ULl, UL2 


GRB071010A 


0.98 (31) 


51.56 


22.4 








s!pULLUL2 


GRB070802 


2.45 (32) 


52.16 


14.7 








S, PUL1,UL2 


GRB090429B 


9.2 (33) 


53.42 


5.8 


GPL 


-0 69+"'" 


46-17+10 72 0.059 7.55 


PUL1,UL2 


GRB081228 


3.4 (34) 


52.43 


3.8 


PL 


1 QQ+0.3I 


0.028 4.44 


PUL1,UL2 


GRB080516 


3.6 (35) 


52.98 


6.8 


PL 




0.039 7.83 


PUL1,UL2 


GRB091221 


< 3.3" 


54.21 


69.0 


PL 


-L62+ii-^^ 


0.043 101.31 


ULl, UL2 


GRB090904B 


< 5.0" 


54.94 


58.2 


PL 


1 cg+8:8i 


0.105 86.40 


UL1,UL2 


GRB 100205 A 


12.0^ 


54.2 


32.76 


PL 




0.008 41.40 


UL2 



Notes. If not mentioned otherwise, the redshift is taken from lGreiner et al.) ( 1201 lah . *** Both the prompt best-fit model and parameters are taken 
from NBll. They are only listed if the E^^a was calculated manually. S: Spectroscopic, P: Photometric, ULl: Upper limit 1, UL2: Up per limit 
2. These upper limits are treated as being the actual redshift value. Value taken from between 1 1 < z < 13.5, iCucchiara et al.ll2010h . 

Referen ces. All referenc es are taken fromlGreiner e t ig(l201 lal Table 2). (l)IWrgani et alJdMol): (2)ICucchiara et al.l ( l2009l): (3) IChomock"etai] 
• 2009bh : (4) Chen et al. (2009); (5) Fvnbo et al. (2009b); (6) Jakobsson et al. (2009); (7) de Ugarte Postigo et al. (2Q09a); (8) iThoene et all 

(12) 
(17) 
(22) 



2009(); (9) Chornock et al. (2009a); de Uaarte Postigo et al. (2010); (10) de Ugarte Postigo et al. (2009b); (11) Cucchiara e Tal.[ (l2008f) 
Berger& Ranch (2008); (13) D'Elia et al. (2008); (14) D'Avanzo et al. (2008); (15) Bergeretal. (2008); (16) Greiner et al. (2009b): 
^Jakobsson et al. (2008b); (18) (Fvnbo et aL .2009a); (19) Thoene et al. (2008b); (20) Perlev et al. (2008); (21) Fvnbo et al., (i2008a) 



Jakobsson et al. (2008d); (23) Jakobsso n et alj (l2008iJ); (24) jV reeswiik et al. (2008); (25) Thoene et al. (2008c); (26) iThoene et ^TflOOSd): 
(27) Malesani et al. (2008); Guidorzi et al.'(l20 09|); (28)|jak obsson et al. ( 20Q8c ); (29)iGreiner et aL (,2009b.); GO),Ledoux et al.l (l2007h;Fox et al.l 
(2008); (31) Prochaska et al. (2007a); (32) Prochaska et a l. (2007b); Elfas dottir et alj ( 120091) ; (33) iTanvij ( 120101) ; (34) lAfonso et al.l ( 120081) ; 
rKruhleretal.i(i2011b) ; (35) iFilgas et al.i (82008) . 



The length of time of GROND observations (covered by 
our sample), AT, is 3.5 years (September 2007 - March 2010) 
and the sohd angle of the BAT detector is = =0.11 

dBartheknv et a l. 2005). Therefore, for GROND, = 
QLaSiUa . QSwift ^ q qjj ^j^^ probabiUty that the LGRB is 
observable to GROND after it has been detected by Swift). 

The LGRB probability, ?7grb, is related to t/bh and rjcow. For 
the first of these parameters, see Eqn. [13] we assume the mass 
above which a BH can form is taken to be Mbh ^ 30Mq, over 
the possible mass range considered of M^in = 0.1 - M^ax - 



100 M(7) ill a Salp eter IMF, following the prescription used 
lYoon etal.l(l2006f) . 



in 



The final parameter to be fixed is the collimation factor, 
T/cou, being derived from th e jet opening angle. The latte r usu- 
ally ranges between 1°-10° (Frail et aIJl200UICenko etai]|2010h 
implying a collimation facto r range of 1.5 X 1 0~"^ - 1.52 x 10"^ 
(utilising rjcou = 1 ^ cos 0jet; iFrail et alJl2()0Th . A median value 
of Ojei - 5°, T/coU = 3.8 X 10"-^ is assumed throughout. 



7 



J. Elliott et al.: The long y-ray burst rate and the correlation with host galaxy properties 
Table 4: Summary of model parameters used in Eqn.fTTI 



Parameter 


Value or Range 


Step size 


Fixed 


'7coll 


3.8 X 10"'' 




y 




0.14 (1) 


... 


y 


riX-my 


0.57 " 


... 


y 


^redshift 


const. 


... 


y 


An 


0.077 


... 


y 


A T 
Ai 


i.z> 




y 


Mbh/Mq 


30.0 




y 


M„,i„/Mo 


0.1 




y 




100 




y 


(Ml, Ma) 


(7, 12). ..(7 + AM, 12 - AM) 


0.05 


n 




0.1 - 1.8 


0.05 


n 


d 


0-2.9 


0.1 


n 


Mother 






n 



Table 5: Median metallicities and masses for different LGRB 
host galaxy programs. 



Survey 


Z 


Sz" 


log,o(M./Mo) 6m' 




(Zq) 


dex 




dex 


Svensson et al. (2010) 


0.54 


0.65 


9.12 


1.84 


Savaalio et al. (2005) 


0.26 


0.12 


9.32 


0.75 


Mannucci et al. (2011) 


0.61 


0.30 


9.33 


0.64 


Rau et al. (2010)* 


0.13 


0.62 






Rriihler et al. (201J^ 






9.76 


0.38 



Notes. Standard deviation. Damped Lyman-a systems at z ~ 2. 



Notes. Product of fraction of GRBs followed by XRT (2) and number 
of GRBs with an X-ray afterglow (3). Calculated based on the sample 
used, i.e., sub-sample size / full-sample size. 



(1) 



iGreiner et al.' ('2011a'); 



References. 

http : //heas arc . nasa . gov/ docs/ swif t/archive/grb_table/stat 
(3) http : //w ww ■ mpe . mpg . de/ - j cg/grbgen . html | 



4.2. Investigated Parameters 

There are four parameters which were varied simultaneously to 
investigate which combination of these best agree with the ob- 
served LGRB rate. The parameters were: galaxy mass lower and 
upper limits (Ali, AI2) respectively, metallicity upper limit {ei), 
the proportionality power (6) and the missing ORB probability, 
'/other, all of which can be found in Eqn.[TT] As before, the param- 
eter ranges investigated can be found in Table |4] and are briefly 
explained below. 

Approximately 80% of LGRB hosts ar e thought to lie withi n 
a mass range of 10'' ''Mq to 1O'' ''M0 ( Savadio et al.ll2009l) . 
To test this hypothesis, a mass range boundary can be applied 
to the CSFH. This is done by moving the values of Mi and 
M2 to a central value, by equal step sizes, i.e., (A1i,Al2) - 
(7, 12).. .(7 + AM, 12 - AM). The initial boundaries are chosen 
to be close to the limits of th at of the GOODS - Field survey and 
are the same as adopted by iBelczvnski et alJ (l2010i) . The step 
size is set to AM = 0.05. 

The LGRB collapsar model requires a metallicity threshold 
to result in a LGRB. Also, LGRB host surveys show a similar 
preference for low-metallicity, as is indicated in Table |5] where 
we give the median metallicities and masses of host galaxies 
from different surveys. Therefore, the metallicity upper limit, 
i.e., the galaxy metallicity below which the galaxy contributes 
to the CSFH, is set to range from 0.1 - 1.8 Zq (the latter value 
being the limit at which the cuts stop affecting the CSFH in the 
model considered for the mass range chosen) with a step size of 
AZ = 0.05 Zq. 

As mentioned previously, it is possible that there are other 
effects not considered by our model or one of the assumptions 
that we made may be inaccurate. This uncertainty is quantified 
by the parameter 6. This black-box approach was carried out by 
[Kistler et al. (2009), who found a value of 6 - 1.5, which could 
be explained by metallicity cuts or an evolving luminosity func- 
tion. As we have included metallicity effects and mass ranges, 
a positive value of 6 would imply that (i) the LGRB LF is non- 
evolving, (ii) the stellar IMF is different to that considered in 
our model or (iii) the progenitor system is incoiTect. To allow 
for a range of possibilities, 6 is varied from to 2.9 in steps of 
Ad^Q.l. 




4 5 6 

Redshift (z) 

Fig. 5: The resulting CSFH model when using different metal- 
licity upper limits. Blue-dashed: = 1-8, green-bold: eL = 1.0 
and red-dash-dotted: fL = 0.1. Each curve utilises (M\,M2) = 



(7, 12) and 6 = 0. The normalisation, 
CSFH from redshift z = - 00. 



P 



Total 



is the integrated 



The final parameter we consider free is //other, obtained from 
Eqn.[TT| As a result of its unknown nature, it will be determined 
from the model fitting that is described in Sect. 14.31 Unlike the 
static parameters described in the Sect. 14.11 the parameters de- 
scribed in this section are redshift and mass dependent. Different 
combinations of these parameters can produce similar outcomes, 
i.e., a degeneracy and each parameter can have an effect on the 
other. For example, when the mass range is brought inwards, 
for each redshift, it will limit also the metallicities, as galaxies 
grow they will produce more metals. The only way to see which 
combination is prefeiTed is to look for the best-fit to the data in 
a systematic approach. The effect of changing each parameter 
individually is depicted in Figs. |5] |6]and [T] 



4.3. Modelling 

The 3D parameter space is to be investigated in a brute-force 
approach, in the following steps: 

1 . The selected samples are first binned in log-luminosity space 
with a bin size of log (i/erg s^'^ = 0.5 for the spectroscopic 
and photometric samples, and 1 .0 for the upper limit 1 and 
upper limit 2 samples (see e.g., Fig.[2]i. 
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10"^ 




Redshift (z) Redshift (z) 



Fig. 6: The resulting CSFH model when using different 
mass ranges. Blue-dashed: (A1i,Al2) = (7,12), green- 
bold: (MuMi) = (7.9, 11.1), red-dash-dotted: (Mi, M2) = 
(8.8, 10.2). Each curve utilises fL = 1-8 and 6 -Q. The normali- 



sation, p 



Total 



is the integrated CSFH from redshift z = - 00. 




Fig. 8: Examples of LGRBR models in comparison to the spec- 
troscopic sample. The blue-solid line is spectroscopic sample, 
the blue-dashed is the model with input parameters 5 - Q, 
eL = 1.8, {M\,M2) - (7, 12) and the red-dashed line is the same 
but with a metallicity limit of = 0.1. All three histograms 
are normalised to their peak value to highlight their main differ- 
ences. 



favoured as the high-z bins contain low counts and can dom- 
inate the slope of the best-fit). 

The data is then compared to the model predictions using 
a least y^-test iy^ = V <-'°p'"'-'--^°''^".d)- \ rj,j^g ^^^^^ ^^^^^ 

assumed to be Poisson distributed. 

This is then repeated for all of the parameters in the 3 dimen- 
sional space described in Table |4] 



An example of the binned data and the corresponding LGRB rate 
models for two metallicity constraints can be seen in Fig [8] 



0123456789 

Redshift (z) 

Fig. 7: The resulting CSFH model when using different power 
dependences. Blue-dashed: 6 - Q, green-bold: S - 1.5, red- 
dash-dotted: 6 - 2.9. Each curve utilises cl = 1.8, (AI1, AI2) - 
(7, 12). The normalisation, pi"'"', is the integrated CSFH from 
redshift z = - 00. 

2. A non-evolving LGRB LF (Eqn.|9]l is then fit to the luminos- 
ity distribution (corresponding to each sample) using a least 
squares algorithm taken from the 5'c/P3Q library. 

3. The chosen sample is then binned in redshift space (bin sizes; 
spectroscopic: 1.74, photometric: 1.87, upper limit 1: 1.6, 
upper limit 2: 2.0) where bin sizes are chosen to ensure a 
limited number of bins (> 80%) have zero counts. 

4. The LGRB number density model (Eqn. fTTT i is then calcu- 
lated utilising a given (Mi , M2), 6, and LGRB LF deter- 
mined in step 2, using the same redshift bin sizes as in the 
previous step, i.e., N(z\,Z2)- 

5. '7othei- is calculated by taking the median value of the ratio 
of the binned data to the expected model data (a median is 

^ |http : //www ■ scipy . org/ 1 



4.4. Summary of Assumptions 

To summarise, the assumptions made throughout this section 
are: 

- LGRBs produce an X-ray afterglow and are coUimated. 

- LGRBs are formed via the single progenitor collapsar model 
from stars above a mass of 30 Mq. 

- A Salpeter IMF between 0.1 - 100 M0, is assumed. 

- LGRBs have a static normal (Gaussian) luminosity function. 

- Galaxies obey a redshift-evolving mass function. 

- Galaxies lie on the mass-metallicity relation. 

- Downsizing is described by an evolving quenching mass 
(mass upper limit). 



5. Results 

Given the methodology outlined in Sect. |4] a LGRB number 
distribution was generated utilising a 3D parameter space con- 
sisting of metallicity limits, mass ranges and missing redshift 
effects. These models were then compared to the different sam- 
ples described in Sect. 13.1! of varying completeness levels, and 
the best-fit results were determined using a least test. 
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(7.00, 12.00) j 



(7.50, 11.50)! 



i (8.00, 11.00)1 



(8.50, 10.50)1 



(8.95, 10.05)' 




(7.00, 12.00), 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.1 




(a) Spectroscopic (S) Sample, d.o.f. = 3. 



(7.50, 11.50)1 



(8.00, 11.00)1 



(8.50, 10.50) 



(8.95, 10.05) 
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(b) Photometric (P) Sample, d.o.f. = 4. 




(7.00, 12.00) 



(7.50, 11.50) 



i (8.00, 11.00) 
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(c) Upper limit 1 (ULl) Sample, d.o.f. = 5. 



(7.00, 12.00) 



(7.50, 11.50) 



: (8.00, 11.00) 



(8.50, 10.50) 
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(d) Upper limit 2 (UL2) Sample, d.o.f. = 5. 



Fig. 9: The best-fit plots for each sample investigated. Top Panels: The contour plot for (5 = 0, white areas denote values much 
larger than the colour scale shown. Only the contour of 5 = is displayed as the fits begin to get worse when 6 starts to increase. 
The black lines denote the maximum to be an acceptable fit to the 1% level. Bottom Panels: The;^^ values for all of the best-fits 
in 6 space, depicting the worst fits for progressively increasing d. 



5.1. Parameter Spaces 

The results of the fits can be seen in Figs. l9all9bll9c] andl9dl Each 
figure depicts the contour plot for metallicity vs mass of the best- 
fit values for 6 = (upper panels). For all of the distributions 
investigated, the best-fit 6 was zero and so only a single contour 
plot is displayed. For clarity to the reader, the best-fit for all 
6 values considered is shown (see lower panels). All the best-fit 
values can be found in Table |6] The zero value of 6 implies that 
there are no missing redshift effects in the modelling that has 
been used (this is discussed further in Sect.|6]l. 

The four samples, S, P, ULl and UL2 show a preference for 
a CSFH-LGRBR connection with no strong global metallicity 
restrictions. The S and P samples show that no mass range lim- 



itations for the host galaxies are preferred. The ULl and UL2 
samples show a preference for the metallicity limits fL = L7 Zq 
and 1.75 Z0 and the mass ranges (Mi, Mi) = (8.0, 11.0) and 
(7.15, 11.85) respectively. These limits are very lax in compar- 
ison to the nomial measured results of < 0.3 Z(?, (see e.g., 
YoonetaL Il2006t I Salvaterra & Chincarinil l2007t ICampisi et al.l 
201 01: IVirdfi et al.ll201 ll) and the range of (M, . Mo) = (9.4,9.6) 
(Savagli o et al.ll200 9). and are interpreted as no limit. The pre- 
dominant difference in our best-fit model compared to previous 
studies in this area is a result of the differences in the peak of the 
redshift number distribution whereby previous samples usually 
peaked at z ~ 4, and required a metallicity cut-off (see Fig. [8]l. 
As we are comparing this to an experimental data set, the way in 
which a sample is selected is very important (this is discussed in 
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more detail in the next section). We note that there is no accept- 
able fit to the S sample to the 1 % level, but is included for com- 
pleteness. Also, there are models present in the contour plot that 
are still acceptable fits at the 1 % level for the following metallic- 
ity ranges and mass boundaries: (All, = (7, 12) -(7.8, 11.2) 
andeL/Zo = 1.65-1.8,(Mi,M2) = (7.95, 11. 05)-(8.05, 10.95) 
and eL/Z© = 1.65 - 1.8, {Mx,M2) = (7, 12)- (8.25, 10.75) and 
fL/Z© - 1 .45- 1 .8, for the P, ULl and UL2 samples respectively 
(1% acceptance boundaries are depicted as bold-black lines in 
Figs|9blEland|9dli. 

5.2. Completeness Levels 

The differences observed between these results and that of pre- 
vious work can be explained by the redshift completeness of 
the samples that have been considered. To generate samples of 
complete redshift, the common approach is to take LGRBs that 
are above a specific luminosity value (luminosity cut). Such a 
method assumes that for all LGRBs above the chosen luminosity 
cut, the probability of measuring the LGRB redshift is equal and 
independent of redshift. A luminosity cut in itself is a perfectly 
acceptable thing to do as one can implement the same choices 
in the model that is used by placing the same luminosity cuts 
in Eqn. |9] However, as it has been mentioned previously, red- 
shift follow-up is not a consistent process. Ground-based pro- 
grams are usually biased towards following potentially interest- 
ing GRBs (i.e., at high-z), and measurements are easier to at- 
tain in the case of the more luminous, quickly detected bursts 
with smaller telescopes. There remains a large gap, both from 
the spectrosc opic desert in the redshift range z = 1 - 2 (see e.g., 
ISteidel et al.ll2004i) . but also in the mid-range redshift z - 2 - 3 
as the LGRB is not visible to the smaller telescopes and the time 
is not always allocated for such bursts on the large telescopes. 
Over the past few years, this has begun to change and all types of 
LGRB are being targeted for follow-up, leading to highly com- 
plete samples, much like the one in this paper For a LGRBR that 
is biased to galaxies with an eL = 0.1, the LGRB number density 
would need to peak at redshift 4 (see Fig.[8]l. Placing luminosity 
cuts on the sample in this paper, it can be seen that the redshift 
peak remains in the range 1-3 (see Fig.fTOll and even very large 
cuts cannot place it in the redshift bin of 4. We also note that 
for the S sample no acceptable fit was found, to the 1 % level, to 
fulfil a null hypothesis, in comparison to the other three samples. 
This is in agreement with the argument presented, as 2/3 photo- 
metric redshifts and 1 upper limit lie in the redshift z ~ 3, which 
are excluded from the spectroscopic sample. 



parameter, //other, a change in the physical parameters t/coU (and 
thus flcou) and t/bh (and thus Mbh) would not change the shape 
of the resulting LGRBR distribution, since the value rjgrb would 
counteract changes in rjcou and t/bh- Assuming that we have in- 
cluded all the LGRB probabilities and thus //other = 1, the result- 
ing changes in ^coii and Mbh can be quantified by: 



M' 



BH 



''coll 



-(^+1) 
'/other 



J Mbh 



(A (m) dm + M"("+" 



COS"' 1 1 



'/other 



(16) 



(17) 



The rescaled values for all the samples investigated can be found 
in Table Q and show two extreme cases: (i) large progenitor 
masses and (ii) small jet opening angles, for which the first prop- 
erty relies on the LGRB collapsar mechanism. The values of (i) 
and (ii) were initially chosen to be Mbh = 30Mq and 0jet - 5° 
respectively (see Sect.|4]i. For a LGRB to occur, we required the 
formation of a BH which we assume forms a progenitor above a 
specific mass, Mbh- Such large masses of the progenitor imply 
that LGRBs form from direct collapse to a BH rather than SN 
fall back under the collapsar model (see e.g.. rHeger et al.l2003h . 
Smaller jet opening angles of ~ 1°, as mentioned previously, 
are possible but would have effects on, for example, the deri ved 
collimated energies of LGRBs (see e.g.. lRacusin et al.ll2009h . 

Despite the fact that we have considered these parameters 
separately, it is also equally valid to tune each value at once to 
give the same final results. To improve upon this, each proba- 
bility parameter must be measured to more precision before the 
CSFH can be used to constrain any single one of them. Finally, 
it is also possible that one of these parameters could also evolve 
with redshift, but as (5 = was shown to be the best-fit in Sect.|5] 
this again is unlikely to be the case. 

6.2. High-z Predictions of the CSFH 

Galaxy detection drops off very quickly with higher redshifts 
due to instrumental limitations making LGRBs complementary 
probes to studying high-z star formation, provided we under- 
stand the relation between the LGRBR and the CSFH. As a re- 
sult of the empirically determined models outlined in Sect.|2]be- 
ing limited to redshifts z < 3 and the possibility that the CSFH 
flattens out at at high redshift, we modify the CSFH m odel to 
a line ar function for redshifts of z > 3 (similar to .Daigne et al.l 
I2006h : 



6. Discussion 

The possible connection between the CSFH and LGRBR was 
investigated utilising a model of the CSFH and a 3D parame- 
ter space of: mass ranges, metallicity limits and proportional- 
ity laws. This was used in combination with a highly complete 
LGRB redshift measured sample. 



6.1. LGRB Probability 

To calculate the LGRB probability parameter, rig^b (Eqn.fT2li. we 
made a priori assumptions about several values (e.g., black hole 
mass and jet collimation). However, the parameters taken from 
the literature that influence the probability of a star forming a 
LGRB have a possible range of acceptable values (see Sect. 14. lb . 
Due to the inclusion in our calculation of //gib of a normalisation 



Piz) 



piz) 



ifz<3 



p (z = 3) - flz if z > 3 



(18) 



where a is a constant to be determined. The parameter a is set to 
vary from the slope of the CSFH at z = 3 and is increased to a flat 
distribution over the range logjo (fl/M0yr"' Mpc^^^ - 0.4-2.5. 
Each parameter fit was then compared to each LGRB sample and 
the best-fit a can be found in Tableland the resulting change in 
the best-fit CSFH can be seen in Fig.fTTI 

All of the samples show no preference for a CSFH flatten- 
ing at high redshift except the upper limit 1 distribution, which 
shows a preference for a CSFH that remains constant with evolv- 
ing redshift. However, the Ax^ between the modified and un- 
modified CSFH is 0.12, and so both are still feasible solutions. 
Overall, the other three samples show no preference for a flatten- 
ing in the CSFH at higher redshifts. This is in contradiction to 
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Table 6: Best fit properties of all the redshift samples. 



Sample 




Ml 


Ml 


S 


'/other 




X'id.oi. 


a 




^modihed * 
/orhRp- 




Zo 


log,o(M/Mo) 


log,o(M/Mo) 




xlO' 






Mq yr"' Mpc"^' 




xlO' 


S 


1.80 


7.00 


12.00 


0.0 


2.5 ± 


1.7 


18.65/3 


0.063 


18.49/3 


2.5 ± 1.7 


P 


1.80 


7.00 


12.00 


0.0 


3.6 ± 


2.6 


4.81/4 


0.050 


6.92/4 


3.6 ± 0.7 


ULl 


1.70 


8.00 


11.00 


0.0 


7.0 ± 


3.2 


12.98/5 


0.013 


12.86/5 


6.3 ± 1.7 


UL2 


1.75 


7.15 


11.85 


0.0 


4.1 ± 


1.9 


2.19/5 


0.032 


4.12/5 


4.1 ±0.5 



Notes. *** Errors are the standard deviation of the first 3 bins. 



Table 7: Redshift and LGRB probability predictions. 



Sample 


5 < z < 6 


6<z<8 


z> 8 


^coU* 




Mbh * 






% 


% 


% 






Mo 




S,v 


3.7 


1.2 


0.0169 


0.98 


+ 0.36 


99.88 


± 36.78 


S" 


2.7 


0.1 


0.0002 


0.98 


+ 0.36 


99.88 


± 36.78 


p« 


3.7 


1.2 


0.0169 


0.82 


±0.30 


99.91 


± 36.76 




5.6 


3.4 


0.0030 


0.82 


±0.30 


99.91 


± 36.76 


ULl" 


4.9 


1.3 


0.0050 


0.58 


±0.21 


99.95 


± 36.77 


ULl" 


7.8 


9.3 


8.0831 


0.60 


±0.22 


99.95 


± 36.77 


UL2'^ 


3.7 


1.2 


0.0147 


0.75 


±0.28 


99.93 


± 36.76 


UL2'^ 


6.9 


6.9 


2.3077 


0.75 


±0.28 


99.93 


± 36.76 



Notes. Normal CSFH model (Eqn. |5j using the best-fit properties given in Table |6] **" Modified CSFH model (Eonll8t using the best-fit 
properties given in Table|6] Uncertainties are given as Icr deviations, as the propagated errors of Eqn.[76land Eqn.ll7lare dominated by l/j/^jij^j 
and are thus underestimated. 



the work bv lKistler etaP (20091), who showed that CSFH would 
be higher than thought, utilising the Swift LGRB sample and a 
single high-z LGRB. As discussed in Sect. 15 .21 this is understood 
as being a result of the extended efforts of m easuring redshifts of 
LGRBs for z > 6. Also, work caiTied out bv lDaigne et al.l (l2006t) 
show that they required a CSFH that increased with redshift. 
However, this would imply higher star formation rates at the 
early stages of the Universe, which are not seen observationally 
or with simulations. Due to this unphysical nature, they propose 
that the LGRB mechanism is evolving with redshift. Again, due 
to the preference of <5 = in our studies, our analysis suggests 
that there is no evolution of the LGRB mechanism required. 

Given the best-fit models the total fraction of LGRBs that ex- 
ist in each redshift range is predicted for each distribution anal- 
ysed (values can be found in Table |7]i. A fraction of ~ 1.2% of 
the LGRB population existing z > 6 is calculated, similar to the 
value of ~ 1% from dCampisi et aLl bOlO). within the error of 
iGreiner et~an (1201 lab w ho predict 5.5 ±2.8% for z > 5 and con- 
sistent with other works (iPerley et al.f20Q9; Fytibo et al. 2009al). 
While this sounds like a small fraction, i t is much larger tha n the 
corresponding fractions for AGN/QSO dWillott et alj|2010h . 

6.3. Summary of Results and Limitations 

It had been previously thought that the LGRB host sample was 
biased to a specific range of masses with low metallicity. We 
find acceptable fits for a mass range of 10^ - IO'^Mq, which 
would suggest that the LGRB rate is a sensitive measure of the 
faint- and massive-end of the mass function of galaxies. We note 
that the proportions of galaxies are determined by the speci- 
fied GMF and the quenching mass, both describe d in Sect. |2l 
The hig h mass range is in good agreement with .Kriihler et al.l 
(1201 la ), who show that by selecting hosts based on high ex- 
tinction more massive galaxies were found. On the contrary, the 
GMF would imply there are more missing low-mass galaxies 
and that LGRBs are selecting star forming regions, with no bi- 
ases. Secondly, the results of Sect. 15.11 also show that there is 



-a 
B 



logio(i„t/erg ) =49 
logi„(L„,/erg s-i)=50 
loSiii(irat/erg s-')=51 
logioiicut/erg ) =52 



012345678 

Redshift (z) 

Fig. 10: LGRB redshift number distribution. Each distribution 
depicts the sample after LGRBs below a specific luminosity limit 
are removed. It can be seen the redshift peak remains in the range 
1-3. 



no strong, ~ 0.1 - 0.3 Z0 metallicity preference for the host 
galaxy at the 95% redshift-complete level. It should be noted 
that this is a global metallicity and does not reflect directly 
on the LGRB model itself, but on the properties of the LGRB 
hosts as a whole. Other types of cuts have been implemented 
where metallicit y dispersio n is also considered (Niino 201 it 
IWolf & Podsiadlowskil2007l) . however, as our model solution re- 
quires no cuts at all, this would not affect the outcome. Thirdly, 
the 6 parameter chosen to quantify any effects not considered 
within our model favours a value of zero. Such a result implies 
the following: (i) the stellar IMF need not be evolving, (ii) the 
LGRB LF is also redshift independent and (iii) the progenitor 
system is consistent throughout all redshifts considered. 



12 



J. Elliott et al.: The long y-ray burst rate and the correlation with host galaxy properties 



10"^ 















\\ 




CSFR \ 




— S modified \ 




ULl modified \ 



0123455789 

Redshift (z) 

Fig. 11: The bold-blue line is the CSFH model described by 
eqn (|5]l. The green-dashed and the red-dotted lines depict the 
best-fit modified CSFH for the S and ULl distributions respec- 
tively. 



In summary, the results show that LGRBs are occurring in 
any type of galaxy (the galaxy number density obeying an evolv- 
ing GMF of full mass range), with no strong metallicity limits. 
The only requirement is that there is star formation occurring, 
in agreement with simulations that suggest LGRB host galax- 
ies preferentially have high specific star formation rates (e.g., 
ICourtv et al.ll2007tlMannucci et al.ll201 ih and observational evi- 
denc e that LGRBs primar ily occur in regions of high star forma- 
tion d^chteretdMoP)- 

As with any model there are limitations. For low redshifts 
(z < 3) we had used empirically calibrated models, that al- 
lowed the freedom of their parameters to be removed (see e.g., 
Eqns. [T] 12] |3]or|4|i. To make reliable judgements at high red- 
shifts (z > 3), further observations able to study these rela- 
tions at higher redshift are needed or other technique s use d 
such as, Monte Carlo Markov Chai n (e.g.. jVir gili et al. 1201 ll) . 
Principal Component Analysis (e.g.. llshida et al. 20 11|), or sim- 
ulations introduced into our model (e .g., Campisi et al.ll201(il 
Ide Souza et alj20r]lllshida et 31.1201? ). Another important con- 
sideration is the assumptions and simplifications made at the be- 
ginning of the modelling, i.e., the stellar IMF, the LGRB lumi- 
nosity function, the LGRB explosion mechanism and the GMF. 
Many of the parameters kept constant could also be left free. 

The form of the stellar IMF, at present, is still a lively de- 
bated issue and could be modified by changing the slope (a), the 
change over short time scales and also the change over redshift 
(top-heavy; see Sect. 12.61 ). Any of these changes would result in 
a modification of the following models that use an IMF for their 
determination: SFR (Sect.ITT), GMF (Sect.lZ2l) and the LGRBR 
probability (Sect. IZ7b . For example. IWang & Pail fcOl 1) show 
that the CSFH-LGRBR connection requires no constraints if a 
redshift dependent stellar IMF is used. For our given framework 
we require no evolution of the IMF, however, any big changes to 
the form of the stellar IMF would require deeper analysis. 

Secondly, the luminosity function of the LGRB was assumed 
to take the form of a normal (Gaussian function). There are many 
other forms, as mentioned previously, such as Schechter func- 
tions, log-normals and redshift dependent functions that could 
be utilised. Again, each one implies different physics, and will 
naturally influence the final result. 



The minimum mass of a star to form a BH was set to Mbh > 
30 M0, as the collapsar model was assumed for LGRB cre- 
ation. Howeve r, there is a range of possible lower range masses 
(iNomoto et al . 2010) and also secondary mechanisms for BHs 
to generate a LGR B, for example SN-fallback (for a review see 
iFrver et al.ll2007h . There are also different types of progenitor 
systems (other than a WR) that are thought to be possible to gen- 
erate a LGRB, e.g., NS-White Dwarf binaries (T homps o n et al] 
l2009h Helium-Hel ium binaries ( Frver & Heg ed |2005|). Quark 
stars touved et all l2005i) . Be/Oe stars (Martavan et al.' '20101; 
Eldridge et al.j 201 lb. blue- straggl ers (Wooslev & Heger 200q) 
and red-giants dEldridge et al.l201 II) . Different mechanisms have 
been investigated befo re, but not in combinat ion with the param- 
eters of this paper (see lYoung & Frverll2007h . 

Finally, the primary difference of the studies of this paper is 
the redshift completeness of the distribution. The normal way of 
improving completeness levels is to choose LGRBs above a spe- 
cific luminosity to compensate for the limitations of the detector 
at increasing redshift. However, redshift measurements are not 
only dependent on the brightness of the afterglow, and follow-up 
of GRBs is not always consistent for many reasons: LGRB sky 
positioning, weather, satellite location and localisation precision, 
to name but a few. Secondly, until recently, many follow-up pro- 
grams were interested in very low and very high redshift GRBs. 
Such a biasing in combination with the redshift desert, results in 
a deficit, in the z ~ 1 - 5 range, that is not easily removed by 
luminosi ty selections (for an analysis of removing selection cri- 
teria see Coward et al .120081) . These deficits can strongly bias the 
results and give different interpretations of metallicity cuts and 
evolving LGRB LFs (see Sect. 15.2b . Rather than improving com- 
pleteness levels by cuts (e.g., luminosity, time criteria), Bayesian 
inference or other methods, incompleteness should be reduced 
by utilising programs with consistent LGRB follow-up. Such 
programs would require no selection biasing of GRB triggers for 
redshift or host follow-up. This is definitely not an easy aim and 
is not always possible with current telescope over-subscription, 
but many programs have alre ady shown that completeness is an 
important criterion (se e e.g ., iFvnbo et al.ll2008"briGreiner et al.l 
2011a; Kruhler et al. 2Q\M. 

In summary, the CSFH-LGRBR connection has a large num- 
ber of free parameters that are co-dependent in several ways, es- 
pecially if redshift dependencies are incorporated. As a result, 
each parameter cannot be treated independently and should al- 
ways be considered contemporaneously, in a systematic way. 
Since the extrapolation of the individual parameters to higher 
redshifts is very difficult to determine, the most direct way of 
improving the uncertainties is primarily a systematic follow-up 
of all redshifts, that would benefit from a GRB mission which 
has a substantially la rger detection rate than Swift (e.g., GRIPS; 
iGreineret al.ll20Hbh . 



7. Conclusion 

The association of LGRBs with the death of massive stars has 
presented many new opportunities for utilising them as high 
redshift tools. One possibility is using the LGRB rate to trace 
the CSFH to unprecedented redshifts, which is usually challeng- 
ing by conventional methods. To reach such a goal, the manner 
in which the LGRB rates traces the CSFH needs to be known, 
whether it is dependent on galaxy properties or not (e.g., galaxy 
stellar mass, metallicity or stellar initial mass functions). In this 
work we have studied whether any type of dependence is re- 
quired, given new evidence of LGRBs occurring in more mas- 
sive galaxies than had been previously thought in combination 
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2009at 



Greiner et al.l 



with highly complete redshift measured LGRB samples (up to 
95% complete). 

Using a highly complete sample we find best-fit solutions 
that show no preference for a strong metallicity or stellar mass 
constraints. These results imply that the LGRB population has 
no preference on the global properties of their host galaxy other 
than it has active star formation. We also show that our ini- 
tial model does not require additional redshift dependences and, 
therefore, implying that there is no redshift dependence required 
in the LGRB probability or luminosity function. 

The best-fit CSFH models are modified at redshifts > 3 to 
linear functions to investigate the possibility that the CSFH flat- 
tens out. The least of the four modified models shows no pref- 
erence for a flattening of star formation at high redshift. This is 
in contradiction with some other LGRB studies, but the exact 
form of the CSFH at high redshifts is still not settled. We predict 
that above z = 6, ~ 1.2% of all LGR Bs exist, whic h is in agree- 
ment with recent simulations (iCamoisi et < 
calstudies (iPerlev et al]|2009l; Ipynbo et al.1 
1201 lah . 

Our results show that sample biasing and completeness lev- 
els of distributions are of essential importance and cannot al- 
ways be recovered in the standard methods. Such completeness 
can only be achieved by consistent follow-up of LGRBs with no 
preference (or bias) on what LGRB is followed. Once unknown 
(and hard to quantify) biases are introduced, they can have dra- 
matic changes to the interpretation of the data. 
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